Magnetic Semiconductors are Frustrated Ferromagnets 
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Starting from microscopic and symmetry considerations, we derive tlie Hamiltonian describing the 
exchange interaction between the localized Mn spins and the valence band holes in Gai-xMuxAs. 
We find that due to the strong spin-orbit coupling in the valence band, this exchange interaction has 
a rather complex structure and generates a highly anisotropic effective interaction between the Mn 
spins. The corresponding ground state has a finite ferromagnetic magnetization but is intrinsically 
spin-disordered even at zero temperature. 
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Semiconductors resisted for decades the concentrated 
efforts of a large community of researchers who wished to 
turn them into magnets. The benefits would be far reach- 
ing: spin could eventually be used to carry information 
in electronic devices. Unfortunately, any conventional 
method to produce semiconductor-magnet alloys failed 
repeatedly: magnetic materials were insoluble in most 
semiconductors. Not long ago ||], several classes of semi- 
conductor materials finally gave in to a most powerful 
weapon against insolubility: the molecular-beam epitaxy 
(MBE) machine. In Gai^xMuxAs, for example, ferro- 
magnetic transition temperatures as high as Tc ~ 110 K 
have been reported [Q. Here we demonstrate, that the 
victory above the natural tendency of the semiconductor 
not to mix with magnetic materials came with a con- 
siderable price. The hostility of the semiconductor host 
against magnetic dopants did not simply disappear, but 
is ever more present in the highly frustrated magnetic cor- 
relations that remain in these systems down to the lowest 
temperatures. We find, that magnetic semiconductors 
are by no means ordinary ferromagnets, but that the in- 
trinsic frustration pushes these systems into the regime of 
strongly spin disordered ferromagnets which also exhibit 
several features reminiscent of spin glasses. 

We start our analysis by constructing the Hamiltonian 
that describes the exchange interaction between a single 
magnetic ion and the valence band holes. In order to 
be specific, we will discuss in detail the important case of 
Gai-xMuxAs. However, our results are rather general, as 
they qualitatively apply to all those magnetically doped 
II- VI and III-V compounds (and in general to any mate- 
rial with magnetic impurities) which show strong valence 
band spin-orbit interaction. 

GaAs is a direct gap semiconductor with a valence 
band maximum at the center of the first Brillouin zone 
called r point Q . The top of the valence band is formed 
by two two-fold degenerate bands of p-character. These 
two valence bands become degenerate in F. The orbitals 
involved in these bands form here a four-dimensional Fg 
irreducible representation. This four-fold degeneracy is 



due to the strong spin-orbit interaction that couples the 
I = 1 angular momentum of the p-orbitals to the elec- 
tron spin (s = 1/2), thereby producing an effective total 
spin J = Z -I- s 3/2 for the valence holes. Since the 
third p-band with J = ^ — s = 1/2 is separated from the 
two by a large spin-orbit splitting. Ago ~ 340 meV, for 
small hole concentrations it is reasonable to describe the 
valence band in terms of a two-band model [Hj5l : 
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where m is the electron mass and the 7i's are the so- 
called Luttinger parameters ||]. The tensor operators 
Jap and Pai3 {a, /? = x,y, x) are defined as Qap — 
\{QaQp + QpQa) - \5ap^i{Clap}, with Q = p and J, 
referring to the momentum of the electrons and their 
J = 3/2 effective spin. In the above equation the last 
two terms, proportional to 72 and 73, describe the cou- 
pling between the effective spin of the valence hole and its 
orbital motion due to spin-orbit interaction. These terms 
will lift the four-fold degeneracy for non-zero momenta. 

In Gai-xMuxAs the Mn ion is believed to be in the 
Mn^^ configuration, corresponding to a half-filled d-shell 
with a total spin 5* = 5/2 The general form of 

the interaction between the S = 5/2 Mn spin and the 
J = 3/2-pseudospin holes depends on the momentum k 
and k' of the incoming and outgoing holes. However, 
close to the F point this momentum dependence is weak 
and we can approximate the coupling constants by their 
k, k' ^ value at the F point. We now proceed to con- 
struct a microscopic Anderson-type Q crystal field model 
that explicitly takes into account both the local crystal 
field symmetry around the Mn impurity and the strong 
Coulomb and Hund couplings. We find (G. Z. and B. J., 
in preparation) that the dominant part of the interaction 
has the following simple form: 

i/int(R) =GS- J(R) (2) 

with G the exchange coupling, and J(R) the spin den- 
sity of the holes at the position R of the Mn ion. Notice 
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that J in Eq. (||) denotes the total J = 3/2 spin of the 
conduction electrons. The above interaction Hamiltonian 
can also be established using purely symmetry considera- 
tions. In general, more complicated couplings of the from 
~ G'J^a^aaSaa, ~ G" J2a^f3 JapSap, etc. are also al- 
lowed by the local symmetry of the Mn ion. However, 
the magnitude of these couplings turns out to be negligi- 
ble compared to G, due to the relatively weak spin-orbit 
interaction on the Mn ion compared to the crystal field 
splitting of the d- levels. 

Equations and (|^) constitute the fundamental 

equations that describe the intricate interplay between 
the spin-orbit interaction in the valence band and the 
local moments. Although the model above can be re- 
fined to incorporate the third valence band, it already 
captures the most important features necessary to un- 
derstand the properties of the ferromagnetic state and 
related phenomena in magnetic semiconductors. One of 
the key differences between this, and earlier models [p^|j4| 
consists in that we now take into account the spin 3/2 
character of the valence holes. This difference, as we 
show below, turns out to be crucial for revealing the true 
ground state of the system. 

In order to analyze the physical content of Eqs. (||) 
and (|l|) we first determine the effective interaction be- 
tween two Mn ions at positions Ri and R2, follow- 
ing the Ruderman-Kittel-Kasuda-Yoshida (RKKY) pro- 
cedure [Q. While the procedure itself is quite straight- 
forward, an explicit evaluation of a more general type of 
interaction deduced in this way is extremely difficult even 
numerically. Fortunately, for the case of GaAs host one 
can make substantial progress by analytical calculations, 
provided that we rewrite Eq. (fil) as 
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The first two parts of the Hamiltonian are rotationally 
invariant, v = (673 -I- 472)/57i ~ 0.77, and the octupolar 
term can be shown to represent only a small correc- 
tion . Therefore, in leading order we can set (5 = and 
consider only the first two, spherically symmetric terms 



which we will denote as H, 



sp- 



in Eq. 

To diagonalize Hgp we choose the spin quantization 
axis to be in the z-dircction. In this basis the energy of 
plane waves propagating along the z-direction is e^(fc = 
kz) = fc^/2m^t with = rrih = m/^i(l — v) Ki 0.5 m 
and — mi — m/^i{\ -\- v) ^ 0.07 m the heavy and 
light hole masses for /i — ±3/2 and fj. = ±1/2, respec- 
tively. Eigenstates of H^-p propagating in other directions 
can then be constructed by simple rotations. The eigen- 
states of _ffsp are chiral in nature: The spin of the heavy 
holes is quantized along their propagation direction k and 
takes the values J • k = ±3/2. In this new chiral basis 
-ffsp is given by the following simple form: 
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where ^ denotes the creation operator of a hole. In this 
basis the unperturbed ground state |0) consists of two 
Fermi spheres. The sphere with the larger radius contains 
heavy holes and includes about 90 % of the valence band 
holes, while a sphere with the shorter radius is generated 
by the light holes. The price for diagonalizing Hs-p is that 
the exchange coupling in the new basis becomes strongly 
momentum dependent: 

Hi„,(R) = f EE^" 4j"(k,k')cke-(''-''')^ . (5) 

k,k' a 

Here J"(k, k') denotes the operator J"(k,k') = 
£'t(k) J"£)(k'), and I?(k) is the spin 3/2 rotation matrix, 
and V the total volume of the sample. It is precisely this 
k-dependence that generates the delicate magnetic prop- 
erties of Gai-xMuxAs. The spherical symmetry of H^p 
implies that the interaction between two impurity spins 
Si and 5*1 at a distance R — |Ri — R2I is given by 



Kp^AR) Sl' • S'l - A'pe,p(i?) Si ■ Si 



(6) 



where S" and denote the spin components parallel 
and perpendicular to Ri — R2. The form of Hcs is some- 
what similar to that of dipolar interactions as it shows 
explicit dependence on the relative position of the Mn 
impurities. Indeed, the interaction between two Mn ions 
far away from each other is in large part mediated by 
holes propagating along the axis R that connects them. 
Since the majority of the holes are heavy and their spin 
is quantized along the propagation direction, it immedi- 
ately follows that the interaction must be different for 
spin components parallel and perpendicular to R. 

It turns out that the structure of the effective inter- 
action Eq. (^) can be calculated analytically, although 
the details are rather technical and will be reported else- 
where . The dominant part of the interaction comes 
from the heavy hole sector, since this has a much larger 
density of states at the Fermi level than the light hole 
band. The heavy hole contribution to Xpar and Kpcrp 
can be expressed as 

-?^par/porp(^) = 27r €f Oh C'par/pcrplfc-F./i-R) , (7) 

where — Ggh is the dimensionless heavy hole exchange 
coupling, gii is the heavy hole density of states at the 
Fermi energy ep, and kp^h denotes the heavy hole Fermi 
momentum. The dimensionless functions Cparp/pcrp(y) 
are clearly different (see Fig. |^), and in the y 
limit are approximately given by Cperp(y) ~ 1/y and 
Cpa.(y) « l/2y. 

In Gai-xMuxAs only a small fraction / of the Mn ions 
gives a hole into the valence band. Although the exact 
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FIG. 1. Main figure: Spatial dependence of tlie two in- 
teraction parameters Cpcrp and Cpar of Eq. (^). Only the 
contribution of the heavy hole sector is shown. The arrow 
indicates the typical value of y for x ~ 0.05 and / ~ 0.25. In- 
set: Two Mn impurities close to each-other tend to be aligned 
•perpendicular to the a^cis connecting them. If we satisfy the 
bonds marked by circles the other bonds with crosses remain 
unsatisfied, resulting in orientational frustration. 
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FIG. 2. Ground state probability distribution function 
of cos Q = n ■ n. The unit vector n gives the direction 
of the magnetization in the ground state. In the isotropic 
approximation, when A'par = ifpcrp = (2/srpar + A'porp)/3, 
the ground state is almost fully polarized, and thus 
P{cos6) = 5{cos9 — 1). Inset: Magnetization as a function of 
temperature for these two cases. Temperature is measured in 
units of 27r ep QhS^ . For the isotropic interaction the magne- 
tization saturates as T ^ 0. 



value of this fraction is not precisely known, latest ex- 
periments suggest that for x=0.05 Mn concentration 
(corresponding to the highest Tc) this fraction is about 
/ = 0.2 - 0.3 (or kF,h « 0.141/1) and a typical Mn-Mn 
distance is approximately dMn-Mn ~ 12A. Thus, for typ- 
ical nearest-neighbor Mn ions ifperp is larger than ifpar 
and ferromagnetic. However, this ferromagnetic inter- 
action is strongly anisotropic as it tries to align nearby 
pair of Mn spins parallel to the axis connecting them (see 
the illustration in Fig. |^). Since it is clearly impossible 
to simultaneously satisfy each pair of spins, this effect 
induces orientational frustration and influences the mag- 
netic properties of Gai-xMuxAs in a fundamental way. 

We investigated the implications of the anisotropy | |T9| ] 
on the magnetic properties of Gai-xMuxAs by per- 
forming classical Monte Carlo (MC) simulations using 
the effective interaction of Eq. (^. In the simulations 
the Mn spins were replaced by classical angular vari- 
ables, S ^ 5 Jl. The Mn ions were randomly dis- 
tributed on a X X face-centered cubic lattice 
with a probability x — 0.05. N = L/a \s the linear 
extension of the lattice in units of the conventional lat- 
tice constant a — 5. 65 A. To take into account the fi- 
nite mean free path I k, 7A of the valence holes, we 
used an exponential cutoff for the RKKY interaction : 
K / (R) ^ K / iR]e~^/^ 

-'^■par/pcrpV *•/ par/perm 

In the inset of Fig. g we show the magnetization 
M = as a function of temperature. We find that 

a spontaneous magnetization develops at low tempera- 
tures |l^ ] . The transition between the paramagnetic and 
magnetic phase takes place rather smoothly, and then in- 



creases approximately linearly with decreasing tempera- 
ture. Both properties agree qualitatively with the exper- 
iments p5[, and are characteristic to strongly disordered 
magnets |10|] . The spontaneous magnetization, however, 
tends to a value at T = that is much smaller then that 
of a fully polarized ferromagnet. This reduction is clearly 
due to orientational disorder originating in the anisotropy 
of the interaction, and has nothing to do with possible 
antiferromagnetic couplings due to the RKKY oscilla- 
tions of the Mn-Mn interaction. To demonstrate this, we 
repeated the simulations by replacing the interaction in 
Eq. (|^) by its angular average. As shown in Fig. ||, the 
magnetization in this case saturates to its maximal value 
(normalized to unity), and all the Mn spins are fully po- 
larized. 

More information can be obtained about the ground 
state properties by measuring the distribution of the 
product cos 6 = $7^11, where n is the a unit vector par- 
allel to the ground state magnetization. Without the 
spatial anisotropy structure discussed here, P(cos 6) — 
5 {cos 6 — 1), since the spins are fully aligned. As shown 
in Fig. |2[ in the system with the correct exchange in- 
teraction, the quantity cos 9 has a very broad but asym- 
metric distribution. Depending on the actual value of / 
the distribution has more or less weight in the vicinity 
of cos 6' = 1: For / = 0.1 the Mn spins tend to point 
approximately into the direction of the global magneti- 
zation, however, they deviate at the average by an angle 
~ 30 degrees and the magnetization is reduced consid- 
erably by ~ 20%. The ground state becomes more spin 
disordered for larger carrier fractions: The distribution 
has less weight at cos 9 — \ and many of the spins are 
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aligned antiferromagnetically with respect to the global 
magnetization, thus reducing the magnetization by about 
- 50%. 

The large reduction in the magnetization we find here 
has been observed experimentally (ijjl^ ■ The measured 
in-plane saturation magnetization of GaMnAs is about 50 
% less than the value that would correspond to the known 
Mn concentration. Also, cooling down the sample in a 
relatively weak {B < AT) external field results in a 20- 
40% increase in the T = magnetization. Also, the broad 
linewidth of ferromagnetic resonance data is consistent 
with the substantial intrinsic spin-disorder reported here. 




The results of our simulations are consistent with 
a highly spin disordered ferromagnetic ground state 
with spin glass-like behavior. Indeed, we found many 
metastable, macro scopically different local energy min- 
ima of the Hamiltonian, extremely close to the ground 
state, a characteristic property of spin glass systems. The 
precise nature of the ground state can be determined ex- 
perimentally. One of the typical experimental signatures 
of a spin glass state would be the history dependence 
of the high field magnetization in fields parallel to the 
film, or the difference between field cooled and zero field 
cooled susceptibilities. This is a straightforward experi- 
ment, which can be performed with already existing sam- 
ples and apparatus. However, the outcome would provide 
a highly valuable insight into the true order present in 
these magnetically ordered semiconducting systems. 

The intrinsic spin disorder we described above could 
be the reason for various resistance anomalies. Since this 
intrinsic spin disorder produces a large spin scattering 
contribution to the resistivity, it may provide an expla- 
nation to the anomalous magnetoresistance of GaMnAs 
alloys with smaller Tc at temperatures T <^ T^. It can 
also possibly explain why these matierals exhibit a re- 
sistance peak precisely at the ferromagnetic phase tran- 
sition: Due to the spin disorder present in the ground 
state, the magnetic order parameter which appears at Tc 
also involves 'disordered' fiuctuations corresponding to 
finite momenta. These finite momentum components of 
the soft modes may result in a maximum of the resistance 
at Tc- This should be contrasted to any conventional fer- 
romagnet, which exhibits predominantly long wavelength 
soft modes and, therefore shows no peak in the resistivity 
at Te 0. 

In conclusion, we have computed the effective exchange 
interaction strong spin-orbit interaction in its valence 
band. We found that the resulting exchange interaction 
inevitably shows a highly anisotropic structure, which in 
turn generates a strongly spin disordered ground state. 
We have shown that the qualitative features of the exper- 
imentally measured magnetization can be reproduced by 
a straightforward simulation of a spin system governed 
by the effective exchange interaction we derived. Our 
results suggest that the actual experimental systems are 



disordered ferromagnets with features reminiscent of spin 
glasses, and proposed an experimental test of such glassy 
behavior. 
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